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The shapes of light normal nuclei and A hypernuclei are investigated in the (3,7) deformation 
plane by using a newly developed constrained relativistic mean field (RMF) model. As examples, 
the results of some C, Mg, and Si nuclei are presented and discussed in details. We found that 
for normal nuclei the present RMF calculations and previous Skyrme-Hartree-Fock models predict 
similar trends of the shape evolution with the neutron number increasing. But some quantitative 
aspects from these two approaches, such as the depth of the minimum and the softness in the y 
direction, differ a lot for several nuclei. For A hypernuclei, in most cases, the addition of a A hyperon 
alters slightly the location of the ground state minimum towards the direction of smaller 8 and softer 
y in the potential energy surface E ~ (8,7). There are three exceptions, namely, BO, PO and 
4'Si in which the polarization effect of the additional A is so strong that the shapes of these three 
hypernuclei are drastically different from their corresponding core nuclei. 


PACS numbers: 21.60.Jz; 21.80.+a; 27.20.+n; 27.30.+t 


I. INTRODUCTION 


Since the first observation of hypernuclei in 1953 [1], 
a lot of experimental efforts have focused on the study 
of the spectroscopy of hypernuclei, see, for examples, 
Refs. [2-4] for recent reviews. Due to the additional 
strangeness degree of freedom, a hyperon is free from 
nucleon’s Pauli exclusion principle. Thus it can move 
deep inside the nuclei and may serve as an impurity for 
probing many nuclear properties that are not accessi- 
ble by normal methods. The study of hypernuclei can 
also provide detailed and accurate information about the 
hyperon-hyperon (YY) and hyperon-nucleon (YN) inter- 
actions [5-8] which are important not only for the un- 
derstanding of hyper nuclear structure but also for the 
study of hyper matter and neutron stars [9]. 

As an impurity in normal nuclei, a hyperon may induce 
many effects on the core nucleus, such as the shrinkage 
of the size [10-14], the change of the shape which will 
be discussed later, the modification of its cluster struc- 
ture [15], the occurrence of nucleon and hyperon skin 
or halo [15-17], and the shift of neutron drip line to a 
neutron-rich side [17-19]. 

The shape describes in an intuitive way the spatial 
density distribution of a quantum many-body system. 
Most of known nuclei are non-spherical and many are 
well-deformed as manifested by regular rotational spec- 
tra [20]. The shape-driven effect of valence nucleon(s) 
has been extensively studied in nuclear high-spin states, 
see, for examples, Refs. [21-24]. A well known example 
is that a nucleon occupying a high-7 and low-Q orbital 
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would drive the nucleus to a more prolate shape. It is 
expected that the addition of a hyperon may also result 
in a shape polarization effect. Since the additional hy- 
peron is not restricted by the nucleons’ Pauli exclusion 
principle, it tends to occupy the lowest s orbital, thus 
driving the core nucleus to be more spherical. If it oc- 
cupies a p orbital, a hyperon may enhance the nuclear 
deformation [25]. 


The self-consistent mean field models, either the non- 
relativistic [26-29] or the relativistic ones [30-42], have 
been extensively used to reproduce the available hyper- 
nuclear data and/or make various predictions for hyper- 
nuclei. Up to now most of these studies focus on spherical 
systems. The first self-consistent mean field study of hy- 
pernuclei with an axially deformed Skyrme Hartree-Fock 
(SHF) model was finished by Zhou et al. [43]. It was 
found that the core nuclei and the corresponding hyper- 
nuclei have similar deformations with the same sign [43] 
which means that the shape polarization effect of the A 
hyperon is quite small. A further study within the same 
framework but with a microscopic AN force gives sim- 
ilar conclusions [44]. However, a relativistic mean field 
(RMF) study reveals that although in most cases the 
results are similar to the SHF calculations, there are 
indeed several exceptions, for examples, {°C and 2°Si 
whose shapes change dramatically compared to their cor- 
responding core nuclei [45]. The different results between 
the SHF and RMF calculations are attributed to the dif- 
ferent polarization effect of the additional A in these two 
approaches [46]. In this sense, the experimental informa- 
tion related to shapes of hypernuclei is much desirable 
and would be used as a good benchmark for theories. 

The triaxiality is an important shape degree of free- 
dom in many nuclei. In an atomic nucleus with a sta- 
bly triaxial shape, the spontaneous broken chiral sym- 
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metry occurs [47-49] and the wobbling motion is also ex- 
pected [50, 51]. The triaxial deformation plays important 
roles in y-soft nuclei and in nuclei in the transitional mass 
regions [52-54]. The y deformation also changes consid- 
erably the local minima and the saddle point in the po- 
tential energy surface of heavy nuclei [55, 56]. Recently 
Win et al. have studied hypernuclei by using a SHF 
model with the triaxial degree of freedom included [57]. 
It is found that with an additional A hyperon no signif- 
icant change occurs for the nuclear shapes except that 
the potential energy surface becomes softer in the y di- 
rection. 

So far the relativistic description of hypernuclei is only 
restricted to spherical or axially deformed cases. It is the 
aim of this paper to investigate the triaxial deformations 
of the A hypernuclei and the shape polarization effect 
of the A hyperon in the RMF model, as a comparative 
study with the SHF calculations [57]. 

The paper is organized as follows. In Section II we 
briefly introduce the RMF model used in the hypernu- 
clear studies with emphasis on the extension of the RMF 
model to the inclusion of the hyperon and the triaxial- 
ity. In Section III, we present the calculated results for 
C, Mg, and Si isotopes and the corresponding hypernu- 
clei and discuss the shape polarization effect of the A 
hyperon. Finally a summary is given in Section IV. 


II. TRIAXIALLY DEFORMED RELATIVISTIC 
MEAN FIELD MODEL FOR HYPERNUCLEI 


In the relativistic mean field (RMF) model, the 
hadrons interact with each other via the exchange of ø, 
w, p mesons and the photon. For hypernuclei, the RMF 
Lagrangian density can be written as: 


L = Lo + La, (1) 


where Lo is the standard RMF Lagrangian density de- 
scribing the nucleons and the couplings between nucleons 
and mesons [58-62] and £4 is that for the hyperon: 


Lie Wa (iy" On, -MA 
+ Ég 


= JoAT — Gury" Wp) PA 


Tn, ^ o” Oa, (2) 


where ma is the mass of the A hyperon, goa and gwa are 
the coupling constants of the A hyperon with the scalar 
and vector meson fields, respectively. The last term rep- 
resents the tensor coupling between the A hyperon and 
the w field [63]. Quv is the field tensor of the w field 
defined as Quy = b pwy — O,w,. Couplings to the p me- 
son and the photon vanishes for A hyperons which are 
neutral and isoscalar. 

Under the mean field approximation, the single particle 
Dirac equation for A hyperons reads: 


[a P+ p(ma + Sa) + Va + Tal Yari = ciai, (8) 


with the scalar potential SA = goĘg, the vector potential 
Va = gwaw and the tensor potential: 


Ta = -4M6 (a By. (4) 
A 


The potential energy surface (PES) is obtained by the 
constrained self-consistent calculation, 


Ne 4 5 
= (H) + 2. z% ((Qn) = Hin) ) (5) 


where H is the RMF Hamiltonian, Qn’s are the multi- 
pole operators to be constrained and Ne is the dimension 
of the constraining space. The quadrupole deformation 
parameters 6 and y are calculated from the multipole 
moments of the baryon density distributions: 


p= p r e a eL (6) 


V3(Q22) 


y = arctan õn) (7) 


where (Ô20) and (Q22) are the quadrupole moments: 
(Âa) = f ar o B-r), 
(Âa) = | dr o0) (0? - 92). (8) 


For normal nuclei, the triaxially deformed RMF model 
has been developed based on expanding the nucleon 
Dirac spinor in a three-dimensional harmonic oscillator 
(3DHO) basis [64, 65]. As an alternative approach, in 
the present work, the RMF equations are solved in an 
axially deformed harmonic oscillator (ADHO) basis [66]. 
We have modified the DIZ (or RMFAXIAL) code [66, 67] 
in order to allow the triaxial deformation and to include 
the hyperon. The basis wave functions are solutions of a 
Schrodinger equation with an ADHO potential: 


Ea IM (o?r? r+ wz »)) la) = Eala), (9) 


with r = yx? +y? and 


1 
a maT exp (m0) Xs 10 
Ja) = on. CIRR (P) EEP (im) xx, (10) 
where a = (nz,N,,m,s) are the asymptotic quantum 


numbers and x, is for the spin. They are characterized 
by the basis deformation g. These basis states form an 
orthonormal complete set and can be used to expand any 
spinor wave functions irrespective with their symmetries. 

The projection of the total angular momentum on the 
symmetric z-axis K is not conserved due to the break- 
ing of the axial symmetry. The remaining symmetries 
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are discrete ones such as interchanges of the three axes. 
Namely, the system is invariant under the point group 
Də. 

To describe the potentials and densities we use a dis- 


mesh points are selected so that a Gaussian quadrature 
can be applied in the r and z directions, while an equally 
distributed mesh is used for the azimuthal angle. For 
convenience we use the Fourier expansion of the poten- 


cretized three-dimensional mesh in the space. These tials V and densities p: 


co 


f (z,7,0) = folz,r) Fz cos (2n8) , (11) 


Fat Lies 


with f = V or p. In an axially symmetric case only the n = 0 term survives, the present code then returns back 
to DIZ. Most of the formulas for f» are formally the same as the corresponding ones for fo, for example, the matrix 


element of the potential V between two basis states is: 


Vaa’ = A s’) 


= 0s,8! (oe 


rR? 
where 
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TABLE I. The calculated binding energies of ?°Si against Np, 
the number of major shells for the Fermion basis. The pa- 
rameter set PK1 is used. The deformation is constrained to 
8 = 0.4 and y = 10°, 30°, and 50°, respectively. The basis for 
bosons are truncated up to Ng = 20. The unit for energies is 


MeV. 
Nr \ (8,7) (0.4,10°) (0.4,30°) (0.4,50°) 
8 —201.667 —202.270 —200.803 
10 —201.523 —202.093 —200.633 
12 —201.382 —201.952 —200.497 
14 —201.325 —201.907 —200.468 
16 —201.296 —201.881 —200.453 
18 —201.294 —201.878 —200.453 


The Klein-Gordon equations for mesons are also solved 
by the basis expansion method, while the Coulomb field 
is solved by the Green’s function method. 

In order to get a point on the PES with given deforma- 
tion parameters (3,7), the axial deformation parameter 
of the basis is set to be, 

Bp = b cosy. (14) 
The harmonic oscillator basis are truncated up to Np 
fermion shells and Ng boson shells. The convergence of 
our method is checked for the nucleus ?°Si. Table I shows 
the calculated binding energies using different Np. Be- 
cause the time consumption of the code is less affected 
by Ng, we set Ng = 20 which is big enough. Three typi- 
cal points on the PES are chosen, i.e., (3,7) = (0.4,10°), 
(0.4,30°), (0.4,50°). In each column in Table I the same 
calculation is performed with Np = 8, 10, , 18. The 


Nez Meri, nl. ay, A OK K+2n,0 + ÔK'— K- N A N E 


n=1 


nim, m 


(12) 
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TABLE II. The RMF parameter sets used in the calculations. 


NN channel ma (MeV) Ro Ro Ruma 
PKLYI [68] PKI [69] 1115.6 0.580 0620 —1 
NLSH-A [45] NLSH [70] 1115.6 0.621 0.667 —1 


truncation errors are less than 100 keV for Nf > 12 and 
less than 30 keV for Np > 14 for these three points. 
Furthermore, when discussing the energy differences on 
which we focus in the present work, the truncation errors 
may even be less due to the cancellation. Therefore we 
use Np = 14 in the following calculations. 


II. RESULTS AND DISCUSSIONS 
A. Numerical details 


In this work we adopted for the Lagrangian density (1) 
two parameter sets which are listed in Table II. For con- 
venience we give in Table II the coupling constants by 
three dimensionless quantities defined as Ro = goA/Jo, 
Ro = gwa/gu and Rusa = fwara/gwa- The effective 
interaction PK1-Y1 is newly proposed by fitting to the 
experimental single-A binding energies and A spin-orbit 
splitting [68]. Based on the parameter set PK1 [69] in 
the NN channel, PK1-Y1 can reproduce the binding en- 
ergies of hypernuclei very well. For comparison we have 
also made calculations using the parameter set labeled as 
NLSH-A which is based on the NLSH parameter set for 
the nucleon-meson coupling constants [70] and has been 
used in the axially symmetric RMF calculations for the 
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hypernuclei in Ref. [45]. 
We use a BCS scheme with constant gaps for the pair- 
ing. Following Ref. [71] the pairing gaps are taken as: 


An = 4.8/N! MeV, A, =4.8/Z1/3 MeV. (15) 


In Ref. [72], the fission barriers in actinides and super- 
heavy nuclei are calculated by using different pairing 
schemes and it is found that BCS calculations with con- 
stant pairing gaps do not provide an adequate descrip- 
tion of the fission barriers. The reason is if, e.g., a con- 
stant pairing strength G is used, the resulting pairing 
gap changes considerably with deformation because the 
density of single particle levels around the Fermi surface 
does so. This is a very important conclusion, especially 
for the study of fission barriers in heavy nuclei. For the 
shape evolution and shape polarization effect of the A hy- 
peron in light nuclei we investigate here, the pairing does 
not play such a decisive role. In Ref. [45] results from a 
constant pairing gap and a constant pairing strength G 
are compared for the light hypernuclei in the RMF+BCS 
model. It is shown that most of the results are only 
slightly changed and the main features and conclusions 
remain the same. However, if the potential energy surface 
is rather soft, the minimum or minima in the potential 
energy surface would be different when different treat- 
ments of the pairing correlation are used. In such cases 
the configuration mixing effects must be also included, 
as discussed later. 

The center of mass correction is included either phe- 
nomenologically (for the parameter sets NLSH and 
NLSH-A) or microscopically (for the parameter sets PK1 
and PK1-Y1). Note that neither the tensor force nor 
the center of mass correction has significant influence on 
shapes of the nuclei investigated in this work. For normal 
nuclei and single-A hypernuclei, they only shift the PES’s 
by roughly a few MeV as a whole. But when discussing 
the absolute value of the energies, their contributions are 
certainly not negligible. 

Using this triaxially deformed RMF+BCS method, we 
have calculated the PES’s of even-even C, O, Ne, Mg, Si 
and S isotopes. Then by adding one A hyperon, we have 
investigated in details the changes of the PES’s of the cor- 
responding hypernuclei. Next we take carbon and silicon 
isotopes as examples and examine the shape evolution 
in these two isotopic chains and the shape polarization 
effect of A hyperon. Some results of 2°Mg and 34"Mg are 
also given for making a comparison with ?°Si and 3’Si. 


B. Carbon hypernuclei 


The ground state properties of carbon isotopes have 
been studied extensively with the axially deformed 
RMF [73-75] and SHF models [76] and the triaxially de- 
formed SHF model [77]. Here we generalize the RMF 
calculations by considering two additional degrees of free- 
dom: the triaxial degree of freedom for nuclear shape and 


the hyperon degree of freedom. We have performed con- 
strained RMF calculations in the (3,7) plane for even- 
even carbon isotopes from !°C to the drip line nucleus 
22C as well as for the hyper counterparts. Note that the 
study of carbon isotopes with triaxially deformed SHF 
models for normal nuclei [77] and hypernuclei [57] are 
both available, while that with the RMF model is still 
absent. 

The calculated ground state deformation parameters 
are summarized in Table III, together with the root mean 
square radii, binding energies and single A separation 
energies. The single A separation energy is defined as 
the energy difference between a hypernucleus and the 
corresponding core nucleus: 


Be OeECO=£20"C). (16) 


Firstly, the contraction due to the additional A are ob- 
served for all the nuclei studies here, which is a manifes- 
tation of the glue like effect of the hyperon [12]. Secondly, 
the radii of '°C and 4'C are much larger than those of 
the neighbouring nuclei with two more neutrons. This is 
due to that they are close to the proton drip line. We 
mention that the earlier calculations with an axially de- 
formed RMF+BCS model also show a similar trend [73]. 
Finally, because the A hyperon always occupies the low- 
est orbital generated by the mean field, the single A sep- 
aration energy By can be seen as a measure of the depth 
and shape of the potential felt by the hyperon. In gen- 
eral this quantity should increase as the nucleon number 
increasing because the depth of the potential increases. 
This is clearly seen in Table III. 

In Figure 1 we present the calculated PES’s of car- 
bon isotopes, together with those of the corresponding 
carbon hypernuclei with one additional A hyperon. The 
locations of the ground states are denoted by open trian- 
gles. To unify the energy scales, we only show the relative 
energies with respect to the ground state. The contours 
join the points on the PES with the same energy. The 
energy difference between two neighbouring contours is 
0.15 MeV. 

Many nuclei show two minima in both the prolate and 
the oblate sides of the PES’s in axially deformed calcula- 
tions. In the triaxially deformed calculations, at least one 
of them becomes a saddle point. In such cases, the rel- 
ative energy differences among the minimum, the saddle 
point and the spherical configuration are key quantities 
characterizing the PES. We list in Table IV the ener- 
gies of the minimum or the saddle point in the prolate 
Eprolate or oblate sides Eopiate as well as that of the spher- 
ical configuration Espherical, With respect to the ground 
state energy. The results calculated with PK1-Y1 and 
NLSH-A parameter sets are both presented. 

It is convenient for further discussions to define the fol- 
lowing two quantities. One is the deformation energy de- 
fined as the energy difference between the spherical shape 
and the ground state, i.e., Eger = |E'spherical — Egrounal- 
Fae characterizes the driving force to deformation quali- 
tatively. The other is the prolate-oblate energy difference 
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TABLE III. The deformation parameters 6 and y, root mean square radii r and binding energies E of the carbon nuclei and 
hypernuclei calculated with parameter set PK1-Y1. The subscripts n, p, A and tot represent the corresponding quantities for 
neutron, proton, A hyperon, and the whole nucleus, respectively. The single-A separation energies Baca are also presented for 
the hypernuclei. 


Nucleus 6 (deg). m. s. radii (fm) Energies (MeV) 
Bn Bp Êa Brot Yn Yp YA tot Mm Tp TA Ttot —Eca Baca 
ug 0.340 0.255 0.284 0 0 0 2.53 2.87 2.74 56.404 
ae: 0.260 0.216 0.060 0.218 0 0 0 0 2.47 2.79 2.43 2.65 65.652 9.250 
Ro 0.207 0.213 0.210 60.0 60.0 60.0 2.53 2.57 2.55 87.764 
C 0 0 0 0 0 0 0 0 2.46 2.50 2.27 2.47 99.752 11.990 
14) 0 0 0 0 0 0 2.74 2.52 2.65 104.943 
rE 0 0 0 0 0 0 0 0 2.73 2.50 2.37 2.62 117.152 12.209 
g 0.225 0.179 0.211 0 0 0 2.90 2.45 2.74 110.398 
Vo 0.202 0.158 0.068 0.183 0 0 0 0 2.88 2.43 2.40 2.70 123.422 13.024 
180) 0.258 0.272 0.261 35.0 40.6 36.4 3.23 2.57 3.03 114.874 
WE 0.246 0.256 0.106 0.243 34.8 40.5 35.7 36.2 3.20 2.54 2.44 2.97 128.434 13.560 
ng 0.275 0.308 0.282 60.0 60.0 60.0 3.47 2.63 3.24 119.601 
O 0.259 0.291 0.119 0.261 60.0 60.0 60.0 60.0 3.43 2.60 2.47 3.17 133.591 13.990 
2o 0.178 0.268 0.193 60.0 60.0 60.0 3.64 2.65 3.40 120.735 
G 0 0 0 0 0 0 0 0 3.55 2.58 2.44 3.28 135.340 14.605 
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FIG. 1. (Color online) The potential energy surfaces of carbon hypernuclei and the corresponding core nuclei in the (6,7) 
plane, calculated using the PK1-Y1 parameter set. The energies are normalized with respect to the binding energy of the 
absolute minimum. The contours join the points with the same energy. The contour interval is 0.15 MeV. The ground states 
are denoted by open triangles. 
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TABLE IV. The energies of the minimum or the saddle point 
in the prolate or oblate sides, and the energy of the spher- 
ical configuration, with respect to the ground state for car- 
bon nuclei and hypernuclei, calculated with the PK1-Y1 and 
NLSH-A parameter sets. The energies are in MeV. 


PKI-Y1 NLSH-A 
Eprolate oblate E’sphercial Eprolate oblate Esphercial 
ug 0 0.74 0.95 0 0.53 0.70 
AC 0 0.34 0.37 0 0.25 0.29 
ug 2 0 0.26 : 0 0.11 
BC - - 0 - : 0 
my : z 0 z 0 
PA £ : 0 7 0 
wo 0 0.18 0.28 0 0.20 0.34 
O 0 0.12 0.18 0 0.14 0.22 
BC 0.30 0.04 1.95 0.22 0.04 1.86 
PC 0.14 0.05 1.57 0.08 0.14 1.52 
200 998 0 2.88 2.15 0 2.75 
"O 1.83 0 2.28 1.74 0 2.14 
220) : 0 0.54 x 0 0.25 
36 - 5 0 : Š 0 


defined as the energy difference between the minimum 
in the prolate (oblate) side and the saddle point in the 
oblate (prolate) side for an axially deformed nucleus or 
the energy difference between two saddle points in the 
oblate and prolate sides for a triaxially deformed nucleus, 
i.e., Epo = |Eprolate — Eoblate|. Apparently it is only well 
defined if at least one saddle point exists in the prolate 
or oblate sides of the PES. Epo characterizes the softness 
of the PES in the y direction in most cases. That is, the 
larger Epo is, the steeper the PES in the y direction is. 


1. Shape evolution of carbon isotopes 


First let us examine the shape evolution of the car- 
bon isotopes. The isotopic dependence of the deforma- 
tion with triaxiality has been investigated with the SHF 
model in Ref. [77], where the PES’s of even-even carbon 
nuclei in the (3,7) deformation plane are presented and 
discussed. Next one can find that there are some new 
features from the RMF calculations. 

The energy minimum of !°C situated at 8 ~ 0.28 and is 
rather stable against the triaxial distortion. The energy 
difference between the ground state and the saddle point 
on the oblate side Epo is 0.74 MeV, accounts for more 
than 1% of the total binding energy. However, this quan- 
tity is less than 40 keV in the SHF calculations [57, 77], 
indicating that 1°C is rather y-soft from the SHF calcu- 
lations. To examine the parameter dependence of our re- 
sults, we also performed the RMF calculation using the 
NLSH parameter set. The resulting Epo is 0.53 MeV. 
Results for !°C from other RMF parameter sets can also 
be found in Ref. [73] where Epo is always larger than 
0.5 MeV. Thus we conclude that the PES of !°C is ap- 
parently softer in SHF than in RMF in the y degree of 


freedom. 


With two additional neutrons in the sd shell, 12C is 
driven to be oblate deformed. The driving force is so 
weak that the deformation energy Eger is only 0.26 MeV 
and 0.11 MeV from the PK1 and NLSH parameter sets, 
respectively. The SHF calculation gives a Eger as large as 
1 MeV [76, 77]. This nucleus is empirically known to be 
oblate in its ground state from the inelastic scattering ex- 
periments [78-80], which is in consistent with our result. 
It is appropriate to compare the depth of the energy min- 
imum obtained from different mean field models. Most 
of the recent SHF calculations of the carbon isotopes use 
a recipe that the spin-orbit interaction is reduced to 60% 
of its original strength in order to reproduce the oblate 
shape of !?C. This prescription is rather arbitrary and 
it is known that the depth of the energy minimum with 
respect to the spherical configuration is sensitive to the 
spin-orbit interaction (see Fig. 1 in Ref. [76]). In Ref. [46] 
it is shown that using this reduction factor as an ad- 
justable parameter, different shapes of a A hypernucleus 
could be predicted. In contrast with the SHF method, 
the deformation of 1?C is correctly reproduced without 
any adjustment of parameters in the RMF model here 
and in Ref. [73]. It seems that the reduction of the spin- 
orbit interaction introduced in these SHF calculations is 
a little too strong, if the RMF results are reliable. This 
needs to be further explored. 


With the neutron number N = 8 which is magic, 14C 
is predicted to be spherical. When two and four more 
neutrons are added to the sd shell, 16C and !8C turn to 
be prolate and triaxially deformed, respectively. For 18C, 
Epo is 0.18 MeV with the PK1 parameter set and 0.20 
MeV with the NLSH parameter set while it is predicted 
to be 0.61 MeV in the SHF calculation [76]. Contrary 
to the result of 1°C, the energy minimum of '®C is much 
deeper in SHF than in RMF. The softness of the PES of 
16C is also obtained by using three different mean field 
models in Ref. [81]. 18C is the only nucleus with a tri- 
axial deformation in carbon isotopes. There are saddle 
points at both the prolate and the oblate sides. The 
RMF calculation predicts that the oblate one is lower 
than the prolate one by 0.26 MeV, while from SHF calcu- 
lations the prolate one is lower (cf. Fig. 1 here and Fig. 5 
in Ref. [77]). The PES’s of 16C and !8C are rather q- 
soft. The ground state deformations of such nuclei with 
extremely soft PES’s may not be well described in the 
mean field level because the ground state wave functions 
are always a strongly correlated superposition of different 
shapes with nearly the same energies. This suggests that 
further investigations of these nuclei should include be- 
yond mean field effects by using the generator coordinate 
method [82, 83]. 

The shape evolves again to be oblate for 7°C and ??C 
with the latter to be the last bound nucleus within the 
neutron drip line of carbon isotopes. ?°C is strongly de- 
formed and its PES is the steepest one among the carbon 
isotopes investigated here in either £ or y directions. ??C 
is suggested to be a halo nucleus according to the mea- 
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sured large enhancement of the reaction cross section for 
it compared to those for neighboring carbon isotopes [84]. 
We certainly can not reproduce the halo structure for ?2C 
because in the present work a Harmonic Oscillator basis 
is used which is not able to give the large spatial density 
distributions in halo nuclei [85-88]. 

The shape evolution of normal carbon nuclei can be 
roughly explained by examining the shell structure. One 
can see from the Nilsson diagram that 10 and 14 are pro- 
late and oblate magic numbers respectively [89], which 
are responsible for the prolate shape of 16C and oblate 
shape of 72°C. 

In the above discussions we see that the RMF and SHF 
models [77] predict the same trend of the shape evolution 
with the neutron number increasing. But some quanti- 
tative aspects from these two approaches, such as the 
depth of the minimum and the softness in the y degree 
of freedom, differ a lot for some nuclei. 


2. Shape polarization effect of A hyperon in carbon 
hypernuclet 


Next let us discuss the shape evolution in carbon hy- 
pernuclei and the shape polarization effect of the A hy- 
peron. As is indicated in Ref. [46], the influence of a 
hyperon on the PES can be as large as 1 MeV, the shape 
may evolve in a different way compared with normal ones. 

Roughly speaking, the deformations of the carbon hy- 
pernuclei are similar to their corresponding core nuclei, 
with two exceptions 43C and 33C. These two nuclei be- 
come spherical while the corresponding core nuclei !7C 
and ??C are both oblate. The spherical shape of 18C 
has been predicted in an axially deformed RMF cal- 
culation [45]. In our triaxially deformed RMF model 
this is confirmed in the (6,7) deformation plane. As 
we discussed earlier, the deformation energy of !?C is 
quite small. From Fig. 1 one can see that at the oblate 
(y = 60°) edge !?C is very soft in the 6 direction. The 
additional A drives \3C to be spherical. %3C becomes 
spherical but it is also very soft in the 8 direction at the 
oblate edge. The calculation with NLSH-A parameter 
set also predicts spherical ground states for {8C and 23C. 

In the SHF calculations, the PES of 1°C with a prolate 
shape is so soft in the y degree of freedom that one addi- 
tional A drives the shape of 4'C to be oblate [57]. In the 
present work, 4'C is still prolate with a smaller 8 com- 
pared to !°C. The addition of a A hyperon only makes 
the PES of }!C a little softer than that of 1°C. Similar 
situation holds for }’C. For the only triaxially deformed 
carbon nucleus, 18C, it is observed that the A hyperon 
also makes the PES of {fC softer and the y softness in- 
creases more towards the prolate direction. Since ?°C 
is strongly oblate deformed, the addition of a A hyperon 
does not change its shape much. Interestingly, with one A 
added, the PES of {C becomes stiffer around the spher- 
ical minimum compared to its core nucleus 14C. This is 
also due to the spherical-driven effect of the A hyperon. 


TABLE V. The energies of the minimum or the saddle point 
in the prolate or oblate sides, and the energy of the spher- 
ical configuration, with respect to the ground state for sili- 
con nuclei and hypernuclei, calculated with the PK1-Y1 and 
NLSH-A parameter sets. The unit of the energies is MeV. 
The results for ?°Mg and 4"Mg are also presented. 


PKI-Y1 NLSH-A 
Eoblate Ebprolate Espherical Eoblate Eprolate E'spherical 
ZZSi - - 0 - - 0 
339i - - 0 - - 0 
249; 0.63 0 1.17 0.17 0 0.62 
255i 0.40 0 0.77 0.07 0 0.33 
26Si 1.02 0 2.20 0.45 0 1.32 
"Si 0.79 0 1.55 0.41 0 0.87 
289i 0 z 1.04 0 - 0.40 
Vei 0 - 0.10 - - 0 
30S; 0 0.17 0.53 0 0.03 0.32 
si <001 0 0.12 0.05 0 0.14 
s29) 0 0.61 0.66 0 0.49 0.52 
3Si 0 0.41 0.44 0 0.32 0.34 
sgi - - 0 - - 0 
25 Si - - 0 - - 0 
26Mg 0.34 0.04 1.38 0.32 0.01 1.08 
"Me 0.29 0 0.89 0.33 0 0.70 


The softness of the PES of a nucleus in the y direction 
can be measured by Epo. For example, it is 0.74 MeV for 
10C and 0.34 MeV for 4' C, which means that the PES of 
the latter is much softer than the former. However, this 
difference for 4’C is only 0.06 MeV smaller than that for 
the corresponding core nucleus, which indicates a very 
tender change. 

From the above discussion it is seen that the shape 
evolution of carbon isotopes is modified due to the addi- 
tional A hyperon. On one hand, the transition point from 
the deformed shape to the spherical shape is shifted from 
140 in normal nuclei to {3C in hypernuclei. On the other 
hand, an abrupt change from a strongly oblate shape to 
a spherical shape is observed for 23C. Thus the sponta- 
neous symmetry breaking effect in carbon hypernuclei is 
very different from that in normal nuclei. 


C. Silicon hypernuclei 


As another example for the study of the shape evolu- 
tion between two shell closures, we present in Fig. 2 the 
PES’s of silicon isotopes and the corresponding one-A 
hypernuclei in the (6,7) plane. Because in silicon iso- 
topes the proton number is 14 which energetically favors 
the oblate shape, the competition between the neutron 
and proton deformation driving forces may produce var- 
ious types of PES’s. When a A hyperon is added, subtle 
changes are expected. 

We list the energies of the minimum or the saddle point 
at the prolate and oblate sides and the spherical config- 
uration, with respect to the ground state in Table V. As 
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FIG. 2. (Color online) The potential energy surfaces of silicon hypernuclei and the corresponding core nuclei in the (3, y) plane, 
calculated using the PK1-Y1 parameter set. The energies are normalized with respect to the binding energy of the absolute 
minimum. The contours join the points with the same energy. The contour interval is 0.3 MeV. The ground states are denoted 


by open triangles. 
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FIG. 3. (Color online) The potential energy surfaces of 7°Mg 
and 4’Mg in the (8,7) plane calculated using the PK1-Y1 
parameter set. The energies are normalized with respect to 
the binding energy of the absolute minimum. The contours 
join the points with the same energy. The contour interval is 
0.3 MeV. The ground states are denoted by open triangles. 


the mirror nucleus of 7°Si, the results for 7°Mg are also 
presented. For silicon isotopes the parameter dependence 
of the results is a bit larger than that for carbon, but the 
results from these two different parameter sets are still 
in consistent with each other. 

Starting from the spherical magic number nucleus 27Si, 


the shape evolves to be prolate in ?4Si and ?°Si. These 
two nuclei are rather y soft, similar as 16C and !8C. In the 
previous SHF calculations [57] an oblate ground state was 
obtained for 7°Si. Following Ref. [57] we also calculate 
26Mg which is the mirror nucleus of 7°Si. The PES of 
26Me is presented in Fig. 3. The ground state of 7°Mg 
is triaxially deformed. Although the calculated PES’s 
for ?°Mg from the RMF and SHF calculations are both 
flat against the y deformation, our RMF model predicts 
a lower saddle point at the prolate side while the SHF 
method does at the oblate side (cf. Fig. 3 here and Fig. 12 
in Ref. [57]). For hypernuclei %Si and 37Si, the values 
of the deformation parameter ( are a bit smaller than 
those of the corresponding core nuclei, respectively and 
the PES’s become softer along the y direction. 


With two more neutrons in the sd shell, ?°Si turns 
to be oblate with Eger as large as 1.04 MeV. Adding 
one A hyperon results in a shape coexistence. Namely, 
the ground state is still oblate, but the energy of the 
spherical configuration with respect to the ground state 
is lowered to as small as 0.1 MeV. The barrier height be- 
tween the two local minima is only about 0.25 MeV. The 
results discussed above are obtained from the parameter 
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set PK1. Note that in the axially deformed RMF calcula- 
tion in Ref. [45], 39Si is predicted to be spherical with the 
NLSH-A parameter set. When the NLSH-A parameter 
set is used in our calculation, we obtain the same conclu- 
sion as Ref. [45]. As is seen in Table V, for 78Si Euer is 
0.40 MeV, but for 4°Si Eqef vanishes which means that 
49Si is spherical. Therefore the prediction of the shape 
of 4°Si is a bit parameter dependent. 

The PES of 2°Si is almost the softest one among the 
nuclei investigated in the present work. Including an ad- 
ditional A hyperon, the energy of 3'Si is almost irrelevant 
with the deformation if 6 < 0.2. The ground state moves 
from the oblate side to the prolate side. But the shift of 
the minimum does not mean much because the coherent 
superposition of the different shapes should be considered 
for the actual ground state. 3?Si is also oblate and the 
PES is soft near the ground state. The A hyperon softens 
slightly the PES of #3Si. Filling completely the neutron 
sd shell, stable spherical shapes are again obtained for 
34Si and 455i. 


IV. SUMMARY 


We developed a triaxially deformed RMF model for 
hypernuclei. Different from previous RMF calculations 
for normal nuclei, in the present work, the RMF equa- 
tions are solved in an axially deformed harmonic oscil- 
lator (ADHO) basis. The convergence of the calculated 
results against the basis truncation is studied and it is 
shown that a reasonably large ADHO basis is able to pro- 
vide desired accuracy in the triaxial RMF calculations. 

The shapes of C, O, Ne, Mg, Si and S A hypernuclei are 
investigated in the (6, y) deformation plane by using this 
newly developed constrained RMF+BCS model with the 
parameter sets PK1-Y1 and NLSH-A. As examples, the 
results of some C, Mg, and Si A hypernuclei are presented 
and we discussed in details the shape evolution of light 
normal nuclei and hypernuclei and the shape polarization 
effect of the A hyperon. 

It is found that for normal nuclei the present RMF 
model and previous Skyrme-Hartree-Fock models predict 
similar trends of the shape evolution with the neutron 
number increasing. But some quantitative aspects from 
these two approaches, such as the depth of the minimum 
and y softness differ a lot for several nuclei. 

For A hypernuclei, in most cases, the addition of a A 
hyperon alters slightly the location of the ground state 
minimum towards the direction of smaller 8 and softer y 
in the potential energy surface (PES) E ~ (6,7). There 
exist two exceptions in carbon isotopes, namely, 4°C and 
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